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Abstract. The dynamics of nucleon transfer processes in heavy-ion collisions is 
investigated within the dinuclear system model. The production cross sections of nuclei 
in the reactions ‘^^Xe-l-^^^Pb and ^'^®U-|-^'^®Cm are calculated and the calculations 
are in good agreement with the experimental data. The transfer cross sections for 
the ^^Ni-l-^^^Pb reaction are calculated and compared with the experimental data. 
We predict the production cross sections of neutron-rich nuclei ‘®®“‘^^Tb, 

i73-i78jjq^ and isi-issyb based on the reaction It can be seen that the 

production cross sections of the neutron-rich nuclei ‘®®Eu, ‘®®Tb, ‘^'^Ho, and ‘®‘Yb 
are 2.84 /xb, 6.90 /xb, 46.24 /xb, and 53.61 /xb, respectively, which could be synthesized 
in experiment. 
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1 . Introduction 

The production of superheavy nuclei (SHN) have been studied extensively in the past 
three decades. Much progress has been made experimentally Bm and theoretically 
[SHn]. There are mainly two sorts of reaction mechanism to produce SHN, which 
are fusion reactions and multinulcleon transfer reactions. The synthesis of superheavy 
nuclei is motivated with respect to searching the “island of stability”. The heavy nuclei 
produced in fusion reactions are not so neutron-rich. Therefore, it is quite hard to reach 
the center of the “island of stability” through the fusion process. The damped collisions 
between heavy nuclei could be one possible approach. 

The light and medium mass nuclei far from the stability line (neutron- and proton- 
rich) can be produced in fragmentation process or fission of heavy nuclei. The heavy 
neutron-rich nuclei including those located at the “island of stability” are not synthesized 
and studied yet, which could be synthesized by using multinucleon transfer reactions. 

The fundamental mechanisms of the deep inelastic collisions have been investigated 
many years ago [El [13]. In recent years, the production of neutron-rich heavy 
and superheavy nuclei have been discussed by Zagrebaev and Greiner [T4l[T5|. The 
multinucleon transfer processes are also investigated by the time-dependent Hartree- 
Fock (TDHF) approach [T6l[T7] and improved quantum molecular dynamics (ImQMD) 
model Ba. The dinuclear system (DNS) model has been successfully used in 
investigating mechanisms of the synthesis of SHN in fusion reactions [T91 - I23] and the 
transfer processes of the damped collisions [21H26] . 

In this work, the transfer cross sections are calculated within the framework of 
DNS model and compared with the available experimental data. The cross sections of 
transferring nucleons from light to heavy fragments and from heavy to light fragments 
are both investigated. We predict the production cross sections of the unknown neutron- 
rich nuclei through the transfer reaction ^^®Yb-|-^^®U. 

The article is organized as follows. In Sec. [2l we describe the theoretical model. 

The results and discussion are presented in Sec. [HI We summarize the main results in 
Sec. [H 


2. Model description 


The diffusion process is treated along proton and neutron degrees of freedom. The 
distribution probability is obtained by solving a set of master equations in the potential 
energy surface (PES) of the DNS. The time evolution of the mass asymmetry is described 
by the following master equation [23] : 




+ -^1, El, t) Ni, Ei,t)] 


N 
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-Agf{Q{t))P{Z,,N,,E,,t). 


( 1 ) 


Here P{Zi, Ni, Ei,t) is the probability distribution function of the fragment with proton 
number Zi and neutron number Ni with corresponding excitation energy Ei at time t. 
^Zi Ni z[ Ni ~ ^Zi Ni Zi n[ mean transition probability from the channel {Zi, Ni, 

El) to {Z[, iVi, E[) [or (Zi, Ni, Ei) to (Zi, E'i)]- dzi,Ni denotes the microscopic 
dimension corresponding to the macroscopic state (Zi, Ni, Ei). The details of W and 
dzi,Ni are given in Ref. [27]. The sum is taken over all possible proton and neutron 
numbers that fragment Z[ and may take, but only one nucleon transfer is considered 
in the model. Agf describes the quasihssion rate. 

The evolution of the DNS along the relative distance R leads to QF of the DNS. 
The QF rate Aqf{Q{t)) can be treated with the one-dimensional Kramers rate [2^ . 


The QF rate exponentially depends on the QF barrier Bqf{Zi,Ni). The local 
temperature 0(f) is calculated by using Fermi-gas expression 0 = i/EfJa with the 
local excitation energy Ei and the level-density parameter a = AjVl MeV~^. The 
frequency of the inverted harmonic oscillator approximates the potential V 'm. R 
at the top of the quasihssion barrier, and u is the frequency of the harmonic oscillator 
approximating the potential in R around the bottom of the pocket. The F determines 
the friction coefficients. Here, F = 2.8 MeV, = 2.0 MeV, and hw = 3.0 MeV. 

The local excitation energy is dehned as 


El — E, 


Idiss - (D(Zi, Ni, Z2, 7V2, J) - D(Zp, iVp, Zt, Nt, J)) 
( J - M)2 m2 




(3) 


Here, U{Zi, Ni, Z 2 , N 2 ) and U{Zp, Np, Zt, Nt) are the driving potentials of fragment 
Ai and the entrance point of the DNS. E^igg is the excitation energy of the composite 
system, which is converted from the relative kinetic energy loss. M denotes the intrinsic 
angular momentum derived from the dissipation of the relative angular momentum, and 
Qnt is the corresponding moment of inertia. J denotes the initial angular momentum. 
(frei is the relative moment of inertia of the DNS, which is given by (rei = pRm- Em is the 
distance of two colliding nuclei located at bottom of potential pocket. For heavy systems 
with no potential pocket the transfer processes take place at the touching configuration. 

During the diffusion process, the relative kinetic energy will dissipate into the DNS 
system. The local excitation energy is determined by the excitation energy of DNS and 
the PES [23|. The PES of DNS is given as 


D(Zi, Ni, Z 2 , iV2, R, A, /?2, 02, J) = Uld{Zi, Ni) + Uld{Z2, N 2 ) - 

Uld{Z, N) - V,^f{J) + Vcn{Zi, Ni, Z 2 , N 2 , R, A, 02, J). 


(4) 


Uld{Zi, Ni), Uld{Z 2 , N 2 ) and Uld{Z,N) are the binding energies of the fragments Ai 
and compound nucleus A, respectively. The value of 17(Zi, iVi, Z 2 , A^ 2 , -R, A, /92, 0i, 02, J) 
is normalized to the energy of the rotating compound nucleus by Uld{Z, N) -f 
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Figure 1. Driving potentials of the tip-tip (solid line) and side-side (dashed line) 
collisions in the reaction ^^®Xe-|-^°®Pb. 


VcN is the interaction potential of two fragments, which depends on the deformation 
parameter and the orientation of the deformed fragments (^i and 02)- The details of 
VcN are given in Ref. m- The distance R between the center of two fragments is chosen 
to be the value at Rm or the touching conhguration, in which the DNS is assumed to 
be formed. 

The sharing of the excitation energy between primary fragments is assumed to 
be proportional to their masses. The code GEMINI is used to treat the sequential 
statistical evaporation of excited fragments. Subsequent de-excitation cascades of the 
excited fragments via emission of light particles (neutron, proton, and a) and gamma- 
rays competing with hssion process are taken into account, which lead to the hnal mass 
distribution of the reaction products. 

The production cross section of a primary fragment with charge Zi and mass 
number Ai can be calculated as follows: 

Jmax Jmax 

= E <uMiEc.m. J)=Y1 ^cap{Ec.m. J)P{Z,, kli, J, E,, (5) 

J=0 J=0 

where CFcapiEc.m., J) is the partial capture cross section, which dehnes the probability of 
overcoming the Coulomb barrier. Ec.m. is the incident energy in center of mass system. 
The details of calculation of (Tcap{Ec,m.) J) is given in Ref. [23]. Tint is the interaction time 
of the diffusion process at the bottom of the potential pocket, which is dependent on the 
incident energy Ec.m. and the incident angular momentum J. Tint cannot be measured 
directly. In this work. Tint is determined by using the deflection function method [28] . 
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Figure 2. The time evolution of the primary mass distributions of the reaction 
^^®Xe+^°®Pb at Ec.m. = 514 MeV with tp = 10“^^ s. is the interaction time 
calculated by using deflection function method [55]. 


3. Results and discussion 

Figure [Tjshows the driving potentials as a function of the mass number of fragments for 
two cases of tip-tip and side-side orientations in the reaction ^^®Xe-l-^°®Pb. The arrows 
show the positions of the masses of projectile and target in entrance channel. A valley 
can be seen around the nucleus ^osp]^ both orientations, which is because ^ospp jg 
a double magic nucleus. It is also can be seen that the valley is deeper for side-side 
orientation. In this work, we investigate the dissipation process under the potential 
energy surface with tip-tip orientation. 

The time evolution of the primary mass distributions, for the reaction ^^®Xe-t-^°®Pb 
at Ec.m. = 514 MeV, is shown in Fig. O It can be seen that the distribution is around 
the entrance channel with Tint = 5 x 10“^^ s. With the increasing interaction time the 
distributions become wider. Equilibrium is reached at about Tint = 40 x 10“^^ s. The 
distribution with Tint = tdf is also shown, t^f is calculated by using deflection function 
method. It is indicated that the system also can reach equilibrium with Tint = tdf, 
which can be seen from mass asymmetry region. The mass distribution around entrance 
channel for Tint = tdf is higher than that for Tint = 40 x 10“^^ s. This is because 
tdf depends on the angular momentum J. The interaction time decreases with the 
increasing angular momentum. 

In Fig. [3l the calculated production cross sections of isotopes of Po, Rn, and Ra 
in transfer reaction ^^®Xe-|-^°®Pb are shown. The bombarding energy Ec,m. =514 MeV. 
The calculated results are the case of tip-tip collisions, which have the height of the 





Production of heavy neutron-rich nuclei in transfer reactions within the dinuclear system modelQ 



M ass num berof fiagm ents 


Figure 3. The comparison of calculated production cross sections for primary (thin 
lines) and survived (thick lines) isotopes of Po (solid lines), Rn (dashed lines), and 
Ra (dash-dotted lines) with the experimental data [29] in the reaction ^^®Xe-|-^°®Pb at 
Ec.m. = 514 MeV. The solid points correspond to the experimental cross sections for 
^^°Po (square), ^^^Rn (triangle), and ^^"‘Ra (star). 



Figure 4. Transfer cross sections for the ®®Ni-l-^°®Pb reaction at Ec.m. =256.8 MeV. 
The number of transferred protons from ^®Ni to ^*^®Pb (negative number) is indicated. 
The measured cross sections have been reported in Ref. m- 
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Figure 5. Transfer cross sections for the ®®Ni+^°®Pb reaction at Ec.m. =256.8 MeV. 
The number of transferred neutrons from ^°®Pb to ®®Ni (positive number) is indicated. 
The measured cross sections have been reported in Ref. [T^ . 


interaction potential at the tonching confignration with the valne 421 MeV. The larger 
qnsifission rate of snch systems with no potential pocket results in the DNS quickly 
decaying into two fragments. Therefore, the neutron transfer process is governed by 
the driving potential with short interaction time. Within the error bars, the calculated 
results are in good agreement with the experimental data [29] for production of ^^“^Po 
and ^^^Ra, while the cross section of ^^^Rn are underestimated. 

Figure m shows transfer cross sections of the reaction ^^Ni+^'^^Pb at E^.m. =256.8 
MeV. Each panel of Fig. 0] shows cross sections according to the change of proton 
number of projectilelike fragment (PLF) from ®®Ni, as functions of neutron number 
of PLF. The calculated results show a good agreement with the experimental data at 
neutron rich region, while the calculated results underestimate the experimental data 
at neutron lack region. As the transferred proton number decreases, the peak position 
of curve shifts towards larger neutron number, while the peaks of the experimental data 
locates at a neutron number of 30. The reason probably is that we only consider the 
driving potential with tip-tip orientation. 

Figure [5] shows transfer cross sections of the reaction ^®Ni-|-^*^®Pb at E^.m. =256.8 
MeV according to the change of neutron number of PLF from ®®Ni, as functions of proton 
number of PLF. From each panel of Fig. [5] one can see that the calculated results are 
in good agreement with the experimental data except that the proton number is 28 in 
first panel. Therefore, the process of transferring some neutrons from ^osp]^ 
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Figure 6. Cross sections for the formation of isotopes of Einsteinium {Z = 99) (solid 
lines), Fermium {Z = 100) (dashed lines), and Mendelevium {Z = 101) (dash-dotted 
lines) in the reaction ^^®U-|-^^®Cm at Ec.m. = 800 MeV. The thin and thick lines are 
distribution of primary and final fragments, respectively. The experimental data are 
taken from Ref. [30]. The circles denote the unknown nuclei. 


transferring some protons from ®®Ni to can be better described. 

In order to further test the model, the cross sections for the production of heavy 
actinides in damped collisions at Ec.m. = 800 MeV are shown in Fig. |6l 

The interaction time of this system is very short due to the strong Coulomb repulsion. 
One can see a good agreement between theoretical and experimental cross sections [30] 
for production of isotopes of the elements Es, Fm, and Md. The cross sections decrease 
drastically with increasing atomic numbers of the fragments. The nuclei ^®^Md and 
262Md could be synthesized in this reaction with cross sections of about 1.52 nb and 
0.17 nb, respectively. Reaction fragments formed in transfer reactions of heavy ions 
are strongly excited. It can be seen that the survival probabilities of most of primary 
fragments are quite low due to a dominant hssion channel and the cross sections of 
primary fragments are a few hundreds times larger than those of the hnal fragments. 
The incident energy should also satisfy the condition that two colliding heavy nuclei can 
come into contact and have enough interaction time for nucleon transfer. Therefore, the 
optimal incident energy should be found for the largest yield of heavy neutron-rich 
nuclei. 

In Fig. [7|the calculated yields of heavy neutron-rich nuclei in the transfer reaction 
i76Yb-|-238u shown. Figure [7| (a) shows the distributions of isotopes of the element 
Yb formed in the reaction at different incident energies. The height of 

the interaction potential at the touching conhguration is about 565 MeV. The height 
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Figure 7. (a) Production cross sections of isotopes of Yb in the transfer reaction 
at Ec.m. = 570 (dash-dotted line), 600 (solid line), 650 (dashed line), and 
690 MeV (dotted line), (b) Cross sections for the formation of isotopes of the elements 
Ytterbium (solid line), Holmium (dashed line). Terbium (dotted line), and Europium 
(dash-dotted line) in the reaction ^’^®Yb-|-^^®U at Ec.m. = 600 MeV. The circles denote 
the unknown neutron-rich nuclei. 


of curve increases with the increasing incident energy from Ec.m. = 570 to 650 MeV, 
while it decreases from Ec.m. = 650 to 690 MeV. This is because with the increasing 
incident energy, the excitation energy of these fragments increases and thus decreases 
their survival probabilities. It can be seen that cross sections are relatively larger in 
neutron-rich region when Ec.m. = 600 MeV. 

The calculated cross sections for production of Europium {Z = 63), Terbium 
{Z = 65), Holmium [Z = 67), and Ytterbium {Z = 70) isotopes in the damped collision 
i76Yb-|-238u qx Ec.m. = 600 MeV are shown in Fig. [7](b). The circles denote the unknown 
heavy neutron-rich nuclei. For the production of neutron-rich nuclei, the cross sections 
decrease dramatically with the increasing neutron number. We show in Fig. [7] (b) 
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the neutron-rich nuclei with production cross section larger than 1 nb. The calculated 
production cross sections for the neutron-rich nuclei I65,i66,i67,i68j^^ 2.84 /xb, 0.78 /xb, 

19.64 nb, and 1.36 nb, I69,i70,i7i,i72,i73rp]^ g_gg 2_';7g gj 2 44.44 nb, and 
1.00 nb, 173 , 174 , 175 , 176 , 177 ,i 78 hq 4 g _24 ^b, 16.83 /ib, 2.00 /ib, 0.85 /ib, 73.75 nb, and 

18.08 nb, and i8i-i82,i83,i84,i85Yb are 53.61 fvh, 4.67 /ib, 2.85 /ib, 84.70 nb, and 4.29 nb, 
respectively. One can see that the cross sections for synthesis of the heavy neutron-rich 
nuclei ^®^Eu, ^®®Tb, ^^^Ho, and ^®^Yb are 2.84 /ib, 6.90 /ib, 46.24 /ib, and 53.61 /ib, 
respectively, which are quite large for experimental detection. 


4. Conclusions 

The calculated production cross sections of heavy nuclei in transfer reaction ^^®Xe-|-^°®Pb 
are in good agreement with the available experimental data. The transfer cross sections 
of the reaction ^®Ni-|-^°®Pb, according to the change of neutron (proton) number of 
PLF from ®®Ni, as functions of proton (neutron) number of PLF are investigated. For 
the production of heavy neutron-rich isotopes of actinides, the reaction 
is also studied. The calculated transfer cross sections show good agreement with the 
experimental data. It is reasonable to use the DNS model to study the dynamics of 
transfer reactions. By using the transfer reaction ^^®Yb-|-^^®U, the neutron-rich nuclei 
^®^Fu, ^®®Tb, and ^®^Yb could be synthesized with the cross sections of 2.84 /ib, 

6.90 /ib, 46.24 /ib, and 53.61 /ib, respectively. 
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